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This paper gives a complete characterization of the location of resonant orbits in a Kerr spacetime 
for all possible black hole spins and orbital parameter values. A resonant orbit in this work is 
defined as a geodesic for which the longitudinal and radial orbital frequencies are commensurate. 

Our analysis is based on expressing the resonance condition in its most transparent form using 
Carlson’s elliptic integrals, which enable us to provide exact results together with a number of 
concise formulas characterizing the explicit dependence on the system parameters. The locations of 
resonant orbits identify regions where intriguing observable phenomena could occur in astrophysical 
situations where various sources of perturbation act on the binary system . Resonant effects may 
have observable implications for the in-spirals of compact objects into a super-massive black hole. 
During a generic in-spiral the slowly evolving orbital frequencies will pass through a series of low- 
order resonances where the ratio of orbital frequencies is equal to the ratio of two small integers. At 
these locations rapid changes in the orbital parameters could produce a measurable phase shift in 
the emitted gravitational and electromagnetic radiation. Resonant orbits may also capture gas or 
larger objects leading to further observable characteristic electromagnetic emission. According to 
the KAM theorem, low order resonant orbits demarcate the regions where the onset of geodesic chaos 
could occur when the Kerr Hamiltonian is perturbed. Perturbations are induced for example if the 
spacetime of the central object is non-Kerr, if gravity is modified, if the orbiting particle has large 
multipole moments, or if additional masses are nearby. We find that the 1/2 and 2/3 resonances 
occur at approximately 4 and 5.4 Schwarzschild radii (Rs) from the black hole’s event horizon. 

For compact object in-spirals into super-massive black holes (~ 10® Mq) this region lies within 
the sensitivity band of space-based gravitational wave detectors such as eLISA. When interpreted 
within the context of the super-massive black hole at the galactic center, Sgr A’*’, this implies that 
characteristic length scales of 41/ras and 55/ras and timescales of 50mm and 79mm respectively 
should be associated with resonant effects if Sgr A’*’ is non-spinning, while spin decreases these 
values by up to ~ 32% and ~ 28%. These length-scales are potentially resolvable with radio VLBI 
measurements using the Event Horizon Telescope. We find that all low-order resonances are localized 
to the strong field region. In particular, for distances r > 50Ra from the black hole, the order of 
the resonances is sufficiently large that resonant effects of generic perturbations are not expected to 
lead to drastic changes in the dynamics. This fact guarantees the validity of using approximations 
based on averaging to model the orbital trajectory and frequency evolution of a test object in this 
region. Observing orbital motion in the intermediate region 50Rs < r < lOOORs is thus a “sweet 
spot” for systematically extracting the multipole moments of the central object by observing the 
orbit of a pulsar - since the object is close enough to be sensitive to the quadruple moment of the 
central object but far enough away not to be subjected to resonant effects. 

PACS numbers: 98.35.Jk, 98.62.Js, 97.60.Lf, 04.20.Dw 


I. INTRODUCTION 

Super-massive black holes such as Sgr A’*’ at the center 
of our galaxy are at zeroth order mathematically ide¬ 
alized as Kerr black holes. In practice this description 
is not complete due to a plethora of small perturbing ef¬ 
fects which slightly alter the spacetime geometry. In gen¬ 
eral these perturbations are small and well accounted for 
with canonical perturbation theory. In the special case 
that the perturbation excites one of the intrinsic reso¬ 
nant structures of the spacetime’s orbits, the effect may 
be larger than normally expected due to an anomalous 
transfer of energy and angular momentum that occurs 


during such a perturbation. Resonance phenomena are 
ubiquitous in any multi-frequency system. In celestial 
mechanics they strongly influence satellite dynamics and 
ring formation. Examples include the gaps in the aster¬ 
oid belt between Mars and Jupiter [l| and the gaps in the 
rings of Saturn M- Resonances are further intimately 
connected with dynamical chaos [^. 

As radio telescopes increase in sensitivity and collect¬ 
ing area we will be able to resolve length-scales typical of 
resonant phenomena in the spacetime of the black hole at 
the center of our galaxy. The Event Horizon Telescope 
is one such observational tool currently under develop¬ 
ment @ . Space based gravitational wave detectors such 
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as eLISA will be sensitive to in-spiral frequencies tra¬ 
verse the resonant regime and may observe shifts in the 
phasing of the gravitational waves emitted during the 
in-spiral of a compact object as it passes through the 
various resonant bands. X-ray, optical and infrared tele¬ 
scopes do not have the resolving power to image Sgr A* 
directly, but can potentially record flux variations from 
this region that may display timescales characteristic of 
resonant events. 

This paper investigates resonant orbits in the Kerr 
metric expanding on the discussion in Q. The aim is 
to provide a complete characterization of the parameter 
space where resonant orbits occur as a function of black 
hole spin and the orbital parameters. Since geodesic or¬ 
bital motion in Kerr is completely integrable, it is akin to 
geodesic flow on a two-dimensional torus in phase space. 
Generic orbits are ergodic and sample the entire surface 
of the torus after a sufficiently long time. Low order reso¬ 
nant orbits however only trace out a simple, co-dimension 
one, curve on the torus. Some of the features of reso¬ 
nant orbital trajectories are illustrated in [Ml. By the 
Kolmogorov-Arnold-Moser (KAM) theorem which is dis¬ 
cussed in Sec. im low order resonant orbits are most likely 
to exhibit the non smooth anomalous behavior associated 
with a rapid change in the constants of motion and the 
breaking of the resonant torus. Test particles entering a 
low order resonance often display subsequent dynamics 
with a sensitive dependence on initial conditions. 

To date a number of authors have studied resonant ef¬ 
fects in Kerr-like metrics in the context of various forms 
of perturbations. The effect of perturbations originat¬ 
ing from adding a quadruple moment to the Kerr metric 
has been quantified by exploring orbital motion in the 
Manko-Novikov metric [l3l - [r^ . Perturbations from the 
presence of a disk were considered in Ref. [l^, and the 
effects of the small mass’ spin in [Tol - I^ . The features of 
traversing a resonance during an in-spiral, where the per¬ 
turbation arises from the small mass’ gravitational self¬ 
force, have been explored by [H, and the possibility 
of sustained resonance has been considered in Res¬ 
onances involving one of the fundamental frequencies of 
the motion on the torus and the orbit’s rotational fre¬ 
quency were studied in the context of enhanced gravi¬ 
tational recoil [1^ [2^ , and isofrequency orbits were dis¬ 
cussed in (T7j . 

Most of these studies have focused on a particular or¬ 
bital trajectory or a small subset of parameters in a spe¬ 
cific perturbed setting. The idea of this paper is to refrain 
from specializing to a particular perturbation and instead 
provide insights that apply to all types of resonant be¬ 
havior. We will use tools such as the results of the KAM 
theorem that hold true regardless of the source of pertur¬ 
bation. The results obtained here are thus robust in the 
sense that the time and length-scales of resonance effects 
for astrophysical applications are to be associated with 
properties of the underlying Kerr metric and resonance 
location rather than the details of the effect causing the 
perturbation. The aim of this paper is to make the typ¬ 


ical resonance time and length-scales accessible to the 
larger astrophysics community by means of easily evalu¬ 
ated formulas and tabulated results. 

To explore the resonance effects we describe the orbits 
in the Kerr metric using a set of variables adapted to 
the orbital geometry that reduce to the Keplerian 
orbital parameters in the Newtonian limit rather than 
the constants of motion associated with the spacetime’s 
Killing fields. The properties of the Keplerian constants 
will be reviewed in Sec. IIII Al Plotting the location of 
resonances in terms of these variables immediately allows 
us to interpret the result as a physical location in the 
actual spacetime. 

A resonant orbit occurs if the ratio of the characteristic 
radial, and longitudinal, wg, frequencies is a rational 
fraction, ujrjojg = n/m where n, m G N. For a more ex¬ 
tended discussion of the orbital geometry and frequencies 
see Sec. IIII Al Most of the technical aspects of this paper 
deal with how to efficiently examine this expression and 
extract the physics. Closed form analytic expressions for 
the frequencies in terms of elliptic integrals have been 
presented by [2^, which serve as companions to this 
work. Here, however, we opt in Sec. IIVI to take advan¬ 
tage of a more symmetric representation of the elliptic 
functions appearing in the resonance condition and write 
them in terms of Carlson’s integrals [sol - ls^ . This al¬ 
lows us to identify the important parameters in the prob¬ 
lem and exploit the identities associated with Carlson’s 
integrals to manipulate the expressions. In Sec. |V] we 
consider solutions to the resonance condition. We first 
specialize to the weak field limit where we introduce the 
key properties of a “resonant surface” in the parameter 
space. We then give a number of exact analytic solutions 
to the resonance condition that can be used to describe 
resonances in the strong field region near the black hole. 
Finally, several low order (small n + m value) resonant 
surfaces such as the 1/2, 2/3 and 3/4 are evaluated nu¬ 
merically and compared to the analytic results and ap¬ 
proximate formulae. 

The breakdown of integrability around a resonance in 
“almost”-Kerr spacetimes is often quantified by numeri¬ 
cally generating Poincare maps for a fixed energy E and 
angular momentum component Lz ■ Associated with each 
Poincare plot is a rotation curve which characterizes the 
frequency ratio as a function of initial condition given 
a fixed E and Lz- In Sec. IVIII we give a representative 
example of orbital breakdown around the 2/3 resonance 
and analytically compute Kerr’s rotation curve. We fur¬ 
ther provide expressions for finding the E and Lz values 
associated with a particular resonance. 

The exact nature of a perturbed system’s response in 
the region of a resonance depends on the source of per¬ 
turbation. In Sec IVIIIl we heuristically discuss how one 
would estimate the size of a perturbation required to see 
a dramatic change in dynamics. It is important to note 
that the KAM theorem does not guarantee the break¬ 
down of integrability at any particular resonance. It 
merely states that if integrability breaks down it will oc- 
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FIG. 1. The orbits in an integrable system with two degrees 
of freedom can be visualized as trajectories wrapping around 
a two dimensional torus in phase space with characteristic fre¬ 
quencies LOi and UJ 2 , relating to the angular advances in 9i and 
62 - For rational values of a; 2 /u;i = mjn the orbital trajectory 
will trace out a distinct path, wrapping n times around the 
9\ axis and m times about the 92 axis. For irrational values 
of L 1 J 2 /wi on the other hand a trajectory will fill the surface of 
the torus densely. 


cur first at the location of a low order resonant orbit 0. 
Since this is true of all possible sources of perturbation, 
the cumulative effect of many sources of perturbation 
could result in a Saturn ringlike structure (see Fig. [HD 
being established around the black hole. This and other 
potentially observable effects due to resonances are dis¬ 
cussed in Sec. m We focus in particular on the galactic 
center, Sgr A*, as a possible observational realization of 
an extreme mass ratio in-spiral (EMRI). We note which 
detectors will be sensitive to chaotic orbits as well as the 
implications of regions were we can guarantee the absence 
of low-order resonances and in which we expect orbits 
to be approximately integrable. Regions that only con¬ 
tain high order resonances we consider to be the “sweet 
spot” for observationally determining the higher order 
multipole moments of the super massive black hole in 
the Milky way. 


II. KAM THEOREM AND IMPORTANCE OF 
RESONANT ORBITS 

Bound geodesic motion in the Kerr spacetimes is in¬ 
tegrable since the Hamiltonian 'Hr = 
where is the inverse Kerr metric and the test par¬ 
ticle’s four momentum, admits a full set of isolating in¬ 
tegrals. Two of these integrals result from the absence 
of explicit time and azimuthal dependence in the Kerr 


* Integrability could also break down at a homoclinic orbit, e.g. 
the last stable orbit discussed in the Appendix However, this 
is of less observational interest than the resonances because it 
marks the transition to the plunge, where the nature of the mo¬ 
tion changes drastically. 


metric functions, the third is due to the conservation 
of rest mass and the fourth integral is known as the 
Carter constant [s^. Integrability implies that action- 
angle variables can be defined. The phase space is foli¬ 
ated by invariant level surfaces of the actions with the 
compact dimensions of these surfaces diffeomorphic to a 
torus. Geodesic motion in an integrable system is thus 
akin to geodesic flow on a torus. 

To illustrate this idea for the Kerr metric, consider the 
reduced Hamiltonian, which is constructed by replacing 
the conjugate momenta associated with the time and az¬ 
imuthal symmetries by their constant values to obtain 
an integrable two degree of freedom system with an ef¬ 
fective potential [l^ . The main features of geodesic flow 
for such a system are sketched in Fig. [TJ The trajectory 
on the torus is described by two characteristic frequen¬ 
cies, associated with the angles 9i and 02 , labeled wi and 
L 02 ^ which in the Kerr metric correspond to the radial 
and longitudinal motions. For rational values of 0 ^ 2 /wi 
the orbit will sample only a finite region of the torus be¬ 
fore retracing its own path, while for irrational values of 
a; 2 /wi a trajectory will fill the torus densely. Orbits with 
rational frequency ratios involving large integers are very 
similar to irrational ones, it is only those with small inte¬ 
ger ratios that are substantially distinct from the ergodic 
case. 

When describing the astrophysical environment 
around a black hole such as Sgr A* we need to take into 
account a number of corrections to the mathematically 
idealized vacuum Kerr metric. In this case we are inter¬ 
ested in the Hamiltonian 


'H = 'Hk + Mi 


( 1 ) 


where "Hi is the perturbing Hamiltonian and e is a di¬ 
mensionless parameter characterizing the strength of the 
perturbation. "Hi contains information about a possi¬ 
ble accretion disk other sources of matter [s^ 

or dark matter (^. l39l| . structural deviations of the cen¬ 
tral black hole away from the Kerr metric (i.e. Bumpy 
Black hole effects [4tH^ b the influence of modified grav- 
ity l^. 1^ . or effects of the multipoles of the small mass 
13 . 4A l48j|. The exact nature of the perturbation does 
not concern us here. In what follows we simply assume 
these modifications to be small and represent this by con¬ 
sidering the case where e <C 1. 

To quantify the effect of an arbitrary perturbation on 
the orbital motion and to find the regions where the 
impact of the perturbation will be greatest, we make 
use of the Kolmogorov-Arnold-Moser (KAM) theorem 
[ 4 ^ [ 5 ^ . The KAM theorem investigates the stability of 
near-integrable systems and suggests that a torus associ¬ 
ated with a rational ratio of frequencies will be destroyed 
in the presence of perturbations. However, provided that 
the perturbation is small enough, tori for which the ratio 
of associated characteristic frequencies are sufficiently ir¬ 
rational will remain stable and persist, although slightly 
deformed, in the perturbed Hamiltonian [SlL l52j | . More 
specihcally, consider the vector of frequencies a; in the 












4 


unperturbed Hamiltonian and a vector of integers k, and 
let d denote the dimension of these vectors. The condi¬ 
tion for resonance is a; • k = 0, which can generally be 
satisfied to arbitrary accuracy by choosing large integers 
for k. When sufficiently large integers are necessary to 
satisfy the resonance condition the tori will be preserved, 
where the definition of sujficiently is such that Arnold’s 
criterion holds 


|a; • k| > K{c) 



( 2 ) 


We will henceforth call Ok = \ki\ the order of the 

resonance. The factor K{e) in Eq. ([2]) approaches zero 
as the perturbation vanishes, i.e. lime_,.o Ar(e) — ?► 0, but 
its functional form depends on the nature of the pertur¬ 
bation. In a non-integrable Hamiltonian system, when 
e < 1, Eq.@ suggests a hierarchy of resonant orbits 
of increasing order whose stability cannot be guaran¬ 
teed. These are the low-order resonances 1, 1/2, 1/3, 
2/3, 1/4, 1/5, 3/4, 2/5, 1/6. We expect these tori to 
be destroyed first if the Hamiltonian is perturbed, how¬ 
ever, from Eq.® we cannot guarantee their destruction 
either. Changing the Kerr metric’s spin parameter is an 
example of a Hamiltonian perturbation to an integrable 
Hamiltonian for which none of the lower order resonant 
tori are broken. 

The destruction of resonant tori corresponds to the 
physical idea that energy transfer takes place most 
rapidly if the frequency of the driving force coincides with 
multiples of the internal frequencies of the system. Sim¬ 
ilarly, even without a direct input of energy, if a system 
is deformed the modes that could potentially be altered 
most are those whose frequencies are rationally related to 
other modes and which thus have the greatest potential 
to exchange energy and interact among themselves. 

The study of torus destruction is not the subject of 
this paper. We do however give a heuristic discussion on 
how to estimate the size of the perturbation required for 
the onset of strongly chaotic dynamics in Sec. IVHII The 
detailed calculation will differ depending on the charac¬ 
teristics of the perturbation. The main focus in the fol¬ 
lowing sections is to identify the regions in parameter 
and physical space where resonant dynamics are likely 
to occur. If they do occur the KAM theorem limits the 
impact to low-order resonances. 


III. GEODESIC MOTION IN THE KERR 
METRIC 

A. Physically motivated constants of motion 

The orbital motion of a bound trajectory of two bod¬ 
ies in Newtonian gravity is described completely by an 
ellipse restricted to a plane. The manner in which this 
ellipse is traversed is characterized by a single frequency, 



FIG. 2. Top Keplerian orbital parameters. The eccentricity, e, is 
a measure of how elliptic the orbit it is. When e = 0, the orbit is 
circular e = 1 the trajectory becomes parabolic. The semi-latus 
rectum, p, can be defined in terms of the eccentricity and the semi¬ 
major axis of the ellipse as is show in green in the figure. The 
point of closest approach rp is called the periastron, while the most 
outlying point the orbit reaches is the apastron denoted by r^. 
Bottom The orbital trajectory as shown in three dimensions. The 
third orbital parameter, namely the maximum inclination angle 
t = 7r/2 — 6*, is the angle with respect to the black hole’s equatorial 
plane and 6, is the minimum Boyer-Lindquist 9 value attained. 


ojcj). A schematic representation of a typical elliptic orbit 
and the Keplerian variables used to describe it is given in 
Fig. [21 By contrast in the Kerr metric bound orbits are 
not restricted to a plane but are confined to a toroidal re¬ 
gion whose shape is characterized by the constants of mo¬ 
tion, the energy E, the 2 —component of angular momen¬ 
tum Lz and Carter constant Q. For geodesics in Kerr, the 
rotational frequency describing the rotational motion 
in the azimuthal direction is augmented by two libration- 
type frequencies ujp and ojg which characterize motion in 
the radial and longitudinal directions respectively. The 
bottom panel in Fig. |2] gives a schematic representation 
of the origin of the uir and wg frequencies associated with 
the orbit. 

In the subsequent sections we will explore the loca¬ 
tion of the resonances for the utr and ujg frequencies 
in Boyer-Lindquist coordinates. Instead of using the 
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constants of motion {E/fi, Lz!p., Q//^^}Q we will de¬ 
scribe the orbits using variables analogous to the Keple- 
rian variables of classical celestial mechanics, namely the 
eccentricity (e), sine of the maximum inclination angle 
(sin i = cos 0 *) and semi-latus rectum (p) as illustrated 
in Fig. [2j These are defined by writing the periastron 
or point of closest approach to the central object as 
Tp =p/(l-|-e), the apastron or furthest point the tra¬ 
jectory reaches as Va =p/(l — e) and the turning point 
of the longitudinal motion as 0*. The typical frequency of 
oscillations between and Vp is described by uj^, whereas 
the longitudinal oscillations about the equatorial plane, 
— (7r/2 — 0*) < t < tt/2 — 0*, are described by wg. 


B. Equations of motion 


For a test mass in orbit around a Kerr black-hole the 
equations governing the radial and longitudinal motion, 
expressed in Boyer-Lindquist coordinates are 



2 


= 



(3) 


where z = cos( 0 ) and we have chosen to parameterize 
the orbit in terms of a non-affine evolution parameter 
A = / dT/ir"^ + a^cos^0), rather than the proper time 
T, so that the radial and longitudinal equations decouple 
(this in fact just corresponds to working in the extended 
phase space). The radial and longitudinal potentials can 
respectively be expressed as 

R = [{r^ + a^)E - aL^] ^ - A -h (L^ - aEf + Q] , 

(4) 

0 = Q(1 - z^) - [(p2 _ E^)a^{l - z^) + Ll] z\ (5) 


where A = — 2Mr + and a = S/M is the spin per 

unit mass. (We will henceforth use units where M = 1.) 
The R and 0 potentials are quartic polynomials of their 
respective arguments and can equivalently be expressed 
in factored form as 


/32 

R = -y(r - ri)(r - r 2 )(r - r 3 )(r - r 4 ) ( 6 ) 


Q = - zl) 


(7) 


where (dS = (/j,^ — E^)a^. In Eq. ([5]) we label the roots so 
that ri > r 2 > ra > 7-4 and in Eq. © so that z+ > Z-. 
For bound orbits, Eq. m dictates that the r(A) and z(A) 
functions describing the orbital motion oscillate between 
two of the roots of Eqs. m and 0 respectively. The 


t The rest mass of the probe, /r, is introduced here to ensure the 
constants of motion are dimensionless. 


generalized Keplerian variables are dehned in terms of 
the roots of the potential functions as: 

ri = , r 2 = , z_ = cos( 0 *) ( 8 ) 

1 — e 1 -I- e 

When quantifying the resonance behavior in the subse¬ 
quent sections we would like to express the results en¬ 
tirely in terms of the variables {p, e, z_} rather than 
using {E, Lz, Q}. The fact that the roots r 4 , and 
z'^ cannot be viewed as independent functions but rather 
must be interpreted as functions of the set of independent 
variables {p, e, z_} complicates the calculation. By com¬ 
paring Eqs. m and © we can find is given explicitly 
in terms of {E, Lz, Q} as follows 

\iLl + Q + /32) ± ^(L2 + Q + ^2)2 _ 4g^2 


(9) 


Equating the coefficients of r in the two expressions 
for the radial equation, Eqs. © and © allows us to 
obtain the following expressions relating the constants 
{E, Lz, Q} to the roots of the factorization 


E^ _ 2(1 -e^) 

p? 2p -I- (1 — e^) w+ ’ 

2 


2 p (p + 2^7+) — 2 a^ (1 — e^) 

2p+ {1 — e^) w+ 

2 {a? (1 — e^) — p^) G7x 

(2p-I-(1 — e^) tx7+) ’ 

Q _ 2p^wx 

p^ (2p -f (1 — e^) t77+) 


( 10 ) 


where we have set vj+ = + r4^ and Wx = In 

addition the condition 


2aELz 


= a 
Pi 


2 (vJx — 2a^) (1 — e^) — 2pwx 
( 2 p+ (1 - e'^)zu+) 


—p^ {w+ — 2 ) + 4ptx7+ 
( 2 p+ (1 - e2)ti7+) 


( 11 ) 


must also hold. Squaring Eq. m and then substituting 
in the expressions for E^ and Eq. m results in a 
quadratic equation for ■cc7_|_ and w- in terms of p and e: 

w+wx p(l - e^)(p -I- 4 - a^) -|-p^(p - 4) - (1 - 

+ vj% [p + - 1 ] ^ + 2pwx [(p - a?) (p -b - l)] -I- 

a'^ (1 — — 2a^ (l — e^) p(p -b 4) -b (p — 4)^p^ 

-b pvj+ [a^ (1 - e^) (a^ + p) - a^P (p -b 4) - (p - 4)p^] 
-bp^ (a^ -p)^ = 0. 


( 12 ) 


If 7 b 0 we can use Eqs. dH]) and (flUl) to rewrite z^ as 


-Z-c = 


P^Wx 


f4(l-e2)z2 ' 


(13) 
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We will always treat the = 0 or Q = 0 limit of orbits 
restricted to the equatorial plane separately. Using Eqs. 
m, cni) and (TT^ we can further show that 

2a^pw+z‘^ = (l — e^) z'^ — P^) ■ (14) 

Eq. (jl4l) is a linear condition in w-\. and Wy_ which, in 
conjunction with (fT^ . implicitly determines '07+(p, e, Z-) 
and ■07x(p, e,z^) and thus the roots and r 4 in terms 
of {p,e,Z-}. Substituting Eq. (fT4)l into (fT^ eliminates 
07+ and yields a quadratic equation in 'cox. As a result 
a closed form expression can easily be found for vjx and 
subsequently ■07+. We will not give the expressions here 
and continue to work with the implicit quantities -07+ 
and -07 X, substituting their actual values only at the end 
of the calculations. The two solutions that result from 
the quadratic equation can be interpreted as test masses 
that either co-rotate or counter-rotate with respect to the 
spin, a, of the black-hole. Orbits that co-rotate with the 
black-hole has the same sign as a) are called prograde 
and those that counter-rotate {L^ has the opposite sign 
to a) are retrograde orbits. For a given {p, e, z_} the pro¬ 
grade orbit’s angular momentum is higher than that of 
the retrograde orbit. On the other hand prograde orbits 
have lower orbital energy than their retrograde counter¬ 
parts 

One special set of orbital parameters is the case when 
the roots satisfy r^ = r 2 , which corresponds to the in¬ 
nermost stable orbit (ISO) separating stable bound or¬ 
bits from those that plunge into the black-hole. For 
a given eccentricity and longitudinal parameter z- the 
semi-major axis satisfying this condition demarcates the 
smallest value of p at which a stable bound orbit can ex¬ 
ist. We shall explicitly solve for the ISO for all values 
of a, e and z- in Appendix and use it as a compara¬ 
tive benchmark for the location of resonant orbits in the 
subsequent sections. 


IV. THE RESONANCE CONDITION 


In this section we begin to characterize the orbits which 
will exhibit resonant behavior. We are interested in the 
parameter values for which uir and ujg are commensurate. 
Given relatively prime integers m and n we seek the sur¬ 
face in the three dimensional parameter space spanned 
by p, e and z_ where 

mojr = nujg. (15) 


This is equivalent to saying that the time it takes the lon¬ 
gitudinal motion to traverse exactly n times between its 
turning points is equal to the time it takes the radial mo¬ 
tion to traverse m times between its turning points. For 
Kerr geodesics m> n since the radial frequency is always 
the smallest of the three frequencies. This translates into 
the following integral condition 


m 



dz _ n dr 

7r' 


(16) 


Substituting Eqs. dH) and © we find that this is equiv¬ 
alent to the condition 

n dr 

an / . = 

dri v(ri - r)(r - r 2 )(r - r 3 )(r-- r4) 



The subject of the rest of the paper is to characterize 
the solutions to this equation. The strategy is to express 
both the radial and longitudinal integrals in their most 
symmetric form using Carlson’s integrals [3 (tI - [^ . Carl¬ 
son’s integral of the first kind is defined to be 

1 dt 

RF{a,P,j) = - , (18) 

2 Jo \/(t + a){t + j 3 ){t + 7) 

In Appendix |B] we list a number of identities and rapidly 
converging approximation techniques that make Carl¬ 
son’s integrals a valuable analytic tool for characterizing 
the resonant surfaces. Using Eq. (IB2I) of Appendix iBl we 
can rewrite Eq. (II3 as 

ani?F(0, (r2 - r3)(ri - r^), (r2 - r4)(ri - rs)) = 

-mi?+(0, (z_-k z+)^, (z_ - z+)^). (19) 

This expression can be further simplified using the iden¬ 
tity ()B8p to rewrite the right hand side and the fact that 
the equations are homogeneous (IB6p to absorb the con¬ 
stant factor. We shall refer to the resulting equation. 


Rf{ 0, {r2 - r3){ri - r4), (r2 - r4)(ri - r3)) = 

i?_F(0, Ka^(z+— z^), Ka^z+), (20) 


as the resonance condition and explore its properties by 
studying various limiting cases. In this expression we 
have defined the parameter 0 < k < 1 to indicate which 
resonance we are considering. 



( 21 ) 


In the subsequent section we will explore all the qual¬ 
itative features of a resonance by examining an easily 
evaluated approximation to Eq. o for large p. We 
then give a number of formulae valid in the region near 
the blackhole in special cases. 


V. SOLUTIONS TO THE RESONANCE 
CONDITION 

When seeking solutions of Eq. (EUl) it is convenient to 
rewrite it in terms of a rapidly converging series. This 
series allows us to identify the three important parame¬ 
ters in the problem. The first sets the overall scale and 
a rough location of the resonance. The remaining two 
are expansion parameters < I that determine the more 
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subtle structure of the resonance surface. We give ex¬ 
plicit expressions for theses parameters in terms of the 
variables introduced in Sec. IIIIBI Next, we evaluate the 
series in the large p limit to obtain a simple analytic 
model which illustrates the important features of any res¬ 
onance. We then turn to the astrophysically more inter¬ 
esting strong field region where the low-order resonances 
occur and give a number of exact analytic formulas for 
special cases. We conclude the section by numerically 
computing the detailed behavior of the 2/3 resonance 
and compare our analytic results and approximations to 
the numerical solutions. 


A. General series expansion 

The resonance condition (l20l) can be rewritten in the 
form 


less than unity ensuring the convergence of the series in 
Eq. (1^ . The parameters Si/yi and ^ 2 / 2/2 vanish when 
e = 0 and Z- = 0 respectively. The special limiting 
case when both these conditions hold, allows us to find 
an exact analytic result that is valid in all regions of the 
spacetime. We shall examine this special case in Sec. IV Cl 
However before we do so, it is instructive to examine the 
properties of resonances that occur at large p values. The 
features we explore in this limit qualitatively capture the 
characteristics of resonances in general. 

B. Anatomy of a resonance in the weak field limit 

p ^ 00 

In the weak field limit, when p —>■ 00 , the dominant 
terms in the expansion oi zux and 737 + found by solving 
Eas. lfHl) and (fl^ are, 


7?f(0, yi+Si,yi- 61 ) = Rf{, 0, K{y 2 + <52), ^( 2/2 - 52))- 

( 22 ) 


where 

p'^-pirs+ri) ep(r4-r3) 

2/1 = - ^5 -H nri, 0 i = — - 5 —, 

1 — 1 — 

2/2 = y(24-^), <^2 = -^. (23) 

It will be shown below Eq. (1251) that 5i/yi <C 1 for all 
physically interesting parameters. As a result, each side 
of Eq. (1221) can be expanded in 5i <C 2/ii using the rapidly 
converging seriesof Eq. (|B10() . Squaring the resulting 
expansions, moving all the terms containing the small 
parameters 5i/yi to one side, and re-expanding the result, 
we obtain the equation 


^=1+1(4-41 

K 2/2 8 \y{ y^J 

123 /Sf 


9 Sfd ^2 

64 yfy^ 


(24) 


In terms of the variables introduced in Sec. Him the 
three quantities that enter the expansion of the resonance 
condition (1^ are, 


2 2 
= a z_ 


Ia7_|_ — 2 • 


1 -I- 2a^(l — z/.) 


± 40-1 




i + -)+o(-). 


(26) 


Substituting Eq. (l26l) into Eq. (l25l) and simplifying the 
result gives the dominant behavior of the three essential 
parameters in the resonance condition 


^^ 6 12a^l - 

2/2 P^ 

3a^((e^-5)z^+4) , 1 \ 

2p2 ^ ^ \p5/2 J 


(27) 




2 e\ 1 — a^z. 


2^2 


4ae 


- 1 


2/1 

P 

^ p3/2 

^ - 

(1 - e2) z4 

■ + o( 

2/2 

2p2 



— I 


O 


(28) 


yi 2a}z^_ ((1 — e^)'DJx + P^ — pi37_|_) 

2/2 -\)zt+ 2 p^Wx) 

-ep^Jw\ - Awx 
yi (I - e2)ti7x-bp2 
82 (1 — e^) 

2/2 (1 — e2)a^zi — 2 p‘^Wx 


(25) 


The first term yi /?/2 ultimately sets the overall scale of p 
at which a particular resonance occurs. Recall that the 
parameters {z_, e, a, k} € [0,1]. By examining Eq. (1251) 
one can further verify that the 8i/yi terms are always 


The small 8 \/y\ and 82 /y 2 parameters scale with 1/p 
and 1 jp^ respectively indicating that they become almost 
negligible for large p and reaffirming the choice of 8 i as a 
expansion parameter. 

To illustrate the basic properties of the resonance con¬ 
dition we substitute Eqs. (1271) and (l28l) into Eq. (l24l) and 
keep terms up to 0{p~^). The resulting approximate res¬ 
onance condition, 

rk28a^Jp{l — z‘/_) = 2 p^{k — 1) -I- I2(p -|- a^) -I- 3e^K 

—3a^z4 (e^(K — 1) -I- 5 ) , 

(29) 

































FIG. 3. Graphical representation of the approximate resonance condition (I29I I or equivalently II32II plotted for the k = (9/10)^ resonance. 
The eccentricity dependence is shown in the three plots on the left. The top left plot displays the typical arch shape seen for all resonances. 
Here the arch is centered around p* ^ 31.6 reaching a maximum value of 2 :_ = 1 at p = Ppoiar given in Eq. (I33II and indicated by a 
dark line on the plot. The lines where arch intersect the z— = 0 plane are indicated in magenta (right) and cyan (left) and show the 
maximum (retrograde) and minimum (prograde) values p attains for a fixed e. Note that the resonance shape for a fixed spin is very 
weakly dependent on the eccentricity of the orbits. The spin dependence of the approximate resonance condition is shown in the three 
plots on the right for e = 1/10 and e = 9/10. Observe that as a —0 the arch pinches off to a line at p = Ppolar-, Fq. Il33|l . The maximum 
arch width occurs for a maximally spinning black hole a = 1 as predicted by Eq. la- 


is valid for large p values only. However, this weak field 
approximation demonstrates all the qualitative proper¬ 
ties of resonant surfaces and gives a good approximation 
even relatively close to the black-hole. The precise man¬ 
ner in which Eq. (1291) breaks down for low order reso¬ 
nances is numerically explored in Sec. IV El 

To build our intuition of the typical features of res¬ 
onant surfaces and their dependence on the parameters 
a, e,Z-,K and p we analyze Eq. (1291) in detail. Eor quasi¬ 
circular orbits and vanishing black-hole spins (a,e) —> 0, 
the resonances occur at 

6 

p* = p{a = 0, e = 0, z_, k) = - -. (30) 

Eor a given k, this value oip* sets the general mean radius 
in physical space (measured in units of GM/(?) about 
which all the interesting features of a resonance occur. 
This is a robust result that remains an exact analytic 
solution even in the region near the horizon, as we will 
prove in Sec. IV Cl For a fixed integer m, resonances with 
K = [(to — 1)/to]^ correspond to the maximum resonance 
radius given by p = 6 to^/(2to — 1). Resonances with n < 
TO — 1 occur at a radius less than that associated with n = 
TO — 1. The maximum p associated with a denominator 
TO thus scales linearly with to for large values of to. 

In the limiting case a —> 0 the dependence on eccen¬ 


tricity is 



Since e < 1 and p* > 6 we see that the effect of eccentric¬ 
ity on the resonance location is small. We also observe 
that in the case a —>■ 0 the location of the resonance be¬ 
comes independent of z_. 

We now examine the general spinning case. Squaring 
both sides of Eq. (1291) results in a polynomial condition 
that is quartic in p and quadratic in We choose 
to analyze the solution surface by specifying the = 
z^(o,e,p, k) for a fixed k rather than explicitly working 
with the quartic roots associated with p. The appropriate 
expression for z^ is 



(32) 
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This function is depicted in Fig. [3] for a fixed spin param¬ 
eter of a = 9/10 and integer ratio k = (9/10)^. It has the 
shape of a parabolic arch centered around p* Ri 31.6. Fur¬ 
thermore, the qualitative features that will be discussed 
here are characteristic for all resonances. The function 
given in Eq. (I32p has a maximum value of z/. = 1 which 
occurs when Ppoiar = p(a, e, Z- = 1, k) has the value 


Ppoiar 



6 (1 — a^) 

p*2 j 

(33) 


Since {e, a} G [0,1] and p* ^ 8, the maximum only de¬ 
viates by a few percent from the p* value as the spin and 
eccentricity deviate from zero. The analytic value of the 
maximum given by Eq. (I33p is plotted as a dark line in 
Fig. [3l When < 1 there are two possible values of p 
that lie on the resonance sheet for a given eccentricity: 
the resonance for a retrograde orbit p_ > Ppoiar and the 
resonance for a prograde orbit which occurs closer to the 
black hole p+ < Ppoiar ■ The sign in the naming conven¬ 
tion of retrograde and prograde orbits relates to the sign 
of the product of the angular momentum and the spin of 
the black hole {aL^), and not the orbit’s relative position 
with respect to Ppoiar ■ 

As decreases and the resonance surface moves from 
the polar towards the equatorial region, the influence of 
spin becomes increasingly important and the distance p_ 
to p-|_ monotonically increases. The expression for p_ 
and p+ can easily be found in closed form by substituting 
Z- = 0 into Eq. and solving the resulting quartic for 
p. However, since the results are messy and add little to 
the discussion we do not give the general results explicitly 
and merely plot these curves in Eig. [31 To benchmark 
the size of the arch we consider the limit of vanishing 
eccentricity and inclination and obtain 


p = Ppoiar (Eq. (l33l) l when a —>■ 0. This indicates that 
resonances in the nonspinning limit become independent 
of inclination because the longitudinal frequency degen¬ 
erates to the ^-frequency in this case. As the black-hole’s 
spin increases from zero the opening width of the arch be¬ 
tween the pro- and retrograde branches increases until a 
maximum arch-width is attained at a = 1. The result 
is a W'-shaped footprint of the arch in the p — a plane, 
with the 'V profile’s vertex corresponding to a = 0. The 
inclination dependence of the resonance surfaces can sim¬ 
ply be characterized as the monotonic closing off of the 
'Y' profile’s pro-and retrograde branches with increasing 
inclination until they merge into a single line forming the 
arch’s spine at Ppoiar- 

This completes our discussion of resonances in the 
weak field limit. The features described here and the 
'U' V I' transitions are characteristic of all reso¬ 
nances. The actual values of the resonant surface of the 
true resonance condition begin to deviate from our weak 
field model as the black hole is approached. The largest 
deviation occurs in the equatorial limit, where the effect 
of spin is most marked. In the polar regions, the weak 
field resonance condition remains a remarkably accurate 
approximation to the true resonance surface. In Sec. IV El 
we numerically characterize several low order resonances 
and provide a quantitative comparison with the approx¬ 
imate results obtained in this section. 

As we shall see next, in the strong field region it is 
possible to obtain exact analytic results for the 'V equa¬ 
torial footprint for e = 0. Since the resonant surface 
depends very weakly on e this result is a good indicator 
for all resonant behavior. 


C. Exact solutions to the resonance condition in 
special cases 


p±(a,e = 0 ,z- = 0, k) _ 1 / [ 4^ 2a \ 

P* 

(34) 

The maximum span of the arch occurs for a maxi¬ 
mally spinning black hole, a = 1. For lower spin val¬ 
ues a good approximation of the span of the arch is 
(p_ — p+) Ri 4ai/^(l + a^/p*). The lowest order cor¬ 
rection to Eq. (I331) with respect to eccentricity is 
e^(p* — 6)/(4p*^) and is the same for both pro and ret¬ 
rograde orbits. 

Having thus explored the basic features of a resonance 
for a given spin parameter a and observed the weak de¬ 
pendence of these features on eccentricity, we will now 
choose a representative eccentricity and then study the 
spin dependence. The right-hand three panels of Eig. |3| 
show the K = (9/10)^ resonance surface for eccentricity 
values of e = 9/10 and 1/10 as a function of black hole 
spin and p. As predicted by Eq. (1341) the arch-width 
exhibits a strong spin dependence. The arch’s inverted 
'U' profile pinches off to a single column profile at 


In this section we explore easily evaluated exact solu¬ 
tions to the resonance condition of Eq. o that can be 
used to characterize the resonant behavior near the black 
hole. The case we will consider first is the limit of circular 
equatorial orbits, i.e. e —?> 0 and —>■ 0. As remarked in 
Sec. IV Al this case sets the parameters di/pi = < 52 / 2/2 = 0 
in Eq. (I33|) and thus a valid solution to the resonance 
condition is found when 


2/2 


= K. 


(35) 


Note that in evaluating this case we will not be resort¬ 
ing to Eq. (I 25 ]) that was derived using Eq. (US which 
assumed that ^ 0. Instead we return to Eq. (jH]) and 
observe that z_ = 0 if and only if Q = 0. By Eq. m we 
see that (5 = 0 implies Wx =0. The simplified version 
of Eq. ([T0|) is 


2 , M=4p-2^^. 

pP 2p -b VO+' pP 2p + VJ+ 


(36) 
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FIG. 4. The location of orbits with resonant frequencies in the limiting case of e = 0 and 2 :_ = 0 as a function of spin, a and 
semi-latus rectum, p [^. Resonances are labeled at their vertex by the rational ratio n/m. All resonances with m < 7 are 
shown. For a = 0 the prograde and retrograde branches are degenerate at p = p*, as the spin increases the retrograde branch 
leans right (copper tinge) and the prograde branch leans left (blue tinge). In general lower order resonances are colored more 
darkly than their higher order counterparts. Note the accumulation of resonances as the strong field region is approached 


Resonance 

Location GMjc^ 

Spin splitting 

Period CMfc'^ 


Galactic center: Sgr A* 

II 

p* = 6/(1 - k) 

Max[(pi - p+)/p*] 

T = 2-kp*'-^/'^ 

p* [pas] T [min] f [10~^Hz] 

ISCO 

6 

1.33 

92.3 

30.6 32.7 5.10 

1/2 

8 

1.22 

142.1 

40.9 50.3 3.31 

1/3 

27/4 = 6.8 

1.29 

110.2 

34.5 39.0 4.27 

2!^ 

54/5 = 10.8 

1.10 

223.0 

55.2 78.9 2.11 

1/4 

32/5 = 6.4 

1.31 

101.7 

32.7 36.0 4.63 

3/4 

96/7 = 13.7 

1.00 

319.1 

70.1 112.9 1.48 

1/5 

25/4 = 6.3 

1.32 

98.2 

31.9 34.7 4.80 

2/5 

50/7 = 7.1 

1.27 

119.9 

36.5 42.4 3.93 

3/5 

75/8 = 9.4 

1.16 

180.4 

47.9 63.8 2.61 

4/5 

50/3 = 16.7 

0.92 

427.5 

85.1 151.3 1.10 

1/6 

216/35 = 6.2 

1.32 

96.3 

31.5 34.1 4.89 

5/6 

216/11 = 19.6 

0.86 

546.7 

100.3 193.5 0.86 

1/7 

49/8 = 6.1 

1.33 

95.2 

31.3 33.7 4.94 

2/7 

98/15 = 6.5 

1.30 

104.9 

33.4 37.1 4.49 

3/7 

147/20 = 7.4 

1.26 

125.2 

37.5 44.3 3.76 

4/7 

98/11 = 8.9 

1.18 

167.1 

45.5 59.1 2.82 

5/7 

49/4 = 12.3 

1.05 

269.4 

62.6 95.3 1.75 

6/7 

294/13 = 22.6 

0.80 

675.8 

115.5 239.1 0.70 


TABLE I. Time and length-scales associated with low-order resonances depicted in Fig. |4]. This table gives the values for the 
e = 0, a = 0, Z- — 0 vertices seen in Fig. [4] first in dimensionless units and subsequently in physical units for the special case 
of the Galactic center, Sgr A*. Lower order resonances, shown in bold in this table, are most likely to have observationally 
detectable dynamics. 


Substituting Eq. (IMl) into Eq. ([9]) gives an expression for Setting Wx = 0 in Eq. dia results in a quadratic equa- 
z+, tion for which has the following roots. 


2 p{p + 2zu+) 
^- 


2p{a ± 

p(p — 4) — =F '^cLy/p 


(37) 


tZ7_|_ — 


(38) 
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We are now in a position to evaluate Eq. (1351) which For low spin values, the solutions to Eq. dm) admit the 
becomes expansion 


m ^ P(P- tZ7+) ^ P-VO+ 

2/2 a?-zX {p + 2w+) 


(39) 


Inserting the value for from Eq. (155)) leads to 


^ {p-6)p-3a^±8a^ 

p^ + 3a2 T 4av^ ’ ^ ’ 

where the upper (lower) sign corresponds to prograde 
(retrograde) orbits. This result can alternatively be ob¬ 
tained from the frequencies of linear perturbations to cir¬ 
cular equatorial orbits . An equivalent way of ex¬ 

pressing Eq. (1401) is in terms of the quartic polynomial; 

[p{p-P*) - a^{p* - 3)]^ - Aa^p{p* - 2)^ = 0. (41) 

In the above expression we chose to use p* = 6/(1 — k) to 
identify the resonance rather than n itself. This choice 
makes it obvious that in the non-spinning limit p = p* is 
an analytic solution to the resonance condition. 

Eq. (HD) is a key result of this paper because it charac¬ 
terizes the exact 'V profile of all resonances for z_ —0 
as a function of spin. As discussed in Sec. IV Bl on the 
weak field limit, eccentricity has very little effect on the 
resonance surface and inclination merely deforms the 'V 
profile into a line as Z- —>■ 1. This single formula thus 
allows us to characterize all resonant effects of arbitrarily 
spinning black holes. 

To efficiently evaluate Eq. dm it is useful to view it as 
a quadratic polynomial in a? instead of a quartic poly¬ 
nomial in p. Solving for in terms of p and p* leads 
to 


_ P^(.P* - 3) -I-p (p* ^ - 5p* -I- 8) 
(p* - 3)2 

-2p(p*-2)v'p(p*-3)-p*+4 
(P* - 3)2 


(42) 


We use Eq. (1421) to plot the spin dependence of the 'V 
profile for = 0 for several low-order resonances in 
Fig. m 

The maximum splitting of the retrograde and prograde 
branches of the 'V occurs when a = 1 in this case the 
relevant roots of Eq. dm are, 


P±=P*-1t2v9^- (43) 


Pzp = p* ± 

a^(p* - 
±— - 


2a(p* — 2) a? (p*^ — 5p* -I- 8) 


p*2 


(45) 


Table U summarizes the numerical values associated with 
the low-order resonances depicted in Fig U) both in di¬ 
mensionless and physical units for the special case of the 
Galactic center, Sgr A*. Lower order resonances, shown 
in bold in this table, are likely to have observationally 
detectable dynamics. According to the KAM theorem 
these tori are most likely to be disrupted and the ensu¬ 
ing rapid changes in the orbital parameters should have a 
dramatic effect when compared to the systematic smooth 
distortion of perturbation induced effects away from res¬ 
onant orbits. We will discuss this further in Sec. m In 
Sec. IVEI we give a numerical characterization of several 
of the lower order resonances introduced here. However 
before we turn to the numerical solution we analytically 
quantify the effect of eccentricity in greater detail. 


D. Quantifying the effect of eccentricity 


To quantify the effect of eccentricity, we consider the 
limiting case of polar orbits with 


= 1, L, = 0, 


2 _ Q _ 

'+~ I3‘^~ a4(I-e2)' 


(46) 


This choice puts us at the top of the inverted 'U' where 
the effects of spin are minimal. In this case we can solve 
Eqs. m and m for vux and vj-. Since is zero 
Eq. (fl^ reduces to a linear equation. As a result, only 
one solution exists as is expected because the retrograde 
and prograde branches should coincide for polar orbits. 

Here to illustrate the effect of eccentricity we only give 
the results in the a = 0 limit, since the expressions for 
a ^ 0 are unwieldy, so that 

Wx =0, w+ = (47) 

p — 4 

Substituting these values into Eq. (1551) gives expressions 
for the parameters that enter into the series expansion of 
Eq. OMl), 


The maximum opening distance of the 'V profile is then 
(P- - P+) = dy'p* - 2. (44) 

Even though Eq. (|4T|) can readily be solved for p, the 
expression is complicated and it is often difficult to iden¬ 
tify which roots correspond to the retrograde and pro¬ 
grade branches, we thus provide a useful series expansion. 


^ _ P- 6 ^ _ 2e 

2/2 P ’ 2/1 p-6’ 


— = 0. (48) 

2/2 


We now re-expand the series in Eq. (l24l) and substitute 
Eq. (|48)) to obtain an explicit expression for k, 


= (I - 6/p) 1 


3e2 

2 (p- 6)2 


516 ^ 

32(p- 6)4 
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FIG. 5. Numerically computed surfaces for the 2/3 resonance. 
Prograde orbits are shown in blue (left) for spin values (a) 0.99, 
(b) 0.5 and (c) 0.2. Retrograde orbits are shown in red (right) for 
spin values (d) 0.2, (e) 0.5 and (f) 0.99. 


Inverting this series expansion we find that the dominant 
behavior of the resonance’s dependence on eccentricity is, 


* 

p = p 



4(p* — 6) 


e'^(4p* - 17) 
64(p* — 6)3 



The expansion given in Eq. ()49|1 is valid in the strong 
field region. Note that as the resonant surfaces approach 
the ISCO (p = 6) the effects of eccentricity become in¬ 
creasingly important. 


E. Detailed numerical characterization of low order 
resonances 

Having studied the resonances in limiting cases, we 
now provide a full numerical characterization of the prop¬ 
erties of the low-order resonances and a comparison with 
the analytic formulas. We quantify the error in using 
the weak field approximation of Sec. IV Bl for low order 
resonant surfaces and confirm the veracity of the exact 
analytic solutions found for special cases in the strong 
field regime. 

We calculated the location of the resonances numer¬ 
ically by means of two methods. The first method di¬ 
rectly uses the closed-form general expressions for the 
fundamental frequencies presented by Schmidt (see 
for example Eqs. [33] and [34] where the frequencies are 
given in terms of elliptic integrals). These formulas are 
not well defined for the limiting cases e = 0 and Z- = 1 
which require explicit modifications [^ . The integrals 
also become indeterminate when reducing to a = 0 or 
a = 1. A careful numerical treatment in the limiting 
cases (e = 0, Z- = 0), where the numerical values are 
calculated using the procedure presented in Appendix B 
of [ 2 ^, recovers the analytic values predicted by Eqs. (|42|) 
and (|49p to within < 0.01% over the entire range of the 
remaining parameters. The accuracy of the agreement in 
these cases are limited by the numerical procedure, e.g. 


taking a = 0.005 and Z- = 0.99995 instead of 0 and 1 
respectively. 

The second semi-analytic scheme that we implemented 
exploits the analytic development of Sec. (IIII Bl) and (IIVI) 
to reduce the computational cost in the following way. 
The resonance condition Eq. (l2^ is recast in terms of 
the three parameters that appear in the series expansion 
of Eq. (l25l) as 



Ug 


/^ i?H0,l + |,l-|) 
V2/2 i?^(0,l + |L,l-|)- 


(50) 


Eor a given a, e, > 0 we calculate the values of p that 
lie on the resonance surfaces by first finding an expres¬ 
sion for CC7+ in terms of vjx using Eq. m and substi¬ 
tuting this into Eq. (El). Solving the resulting quadratic 
for Wx yields two choices for vaxip), corresponding to 
the prograde and retrograde solutions. Having expressed 
Wx and thus in terms of p we can write the three 
parameters in Eq. (1251) as explicit functions of p. The 
right-hand side of Eq. (ISO)) can be evaluated for a given 
p by repeatedly applying the relation (|B8p which tends 
to make the two arguments in the function 

equal (5 iterations usually suffice to reach machine preci¬ 
sion) and then using (IB7I) to yield the result; alternately 
one can use Eq. (IB12P to express Carlson’s Rp functions 
as elliptic integrals of the first kind. A line search method 
can then be used to find a value p that solves (|50|) for a 
given K, using p = p* as the first guess. This method has 
yielded accurate results for all parameter values except 
in the equatorial case Z- = 0. The cases e = 0, z- = 1, 
0 = 0 and a = 1 do not require special treatment, in con¬ 
trast to the first direct numerical method. The special 
case of z_ = 0 can be treated analytically, as was shown 
in the previous sections. 

In the numerical investigation shown in Figs. [5] and [S] 
we sample the parameter space by choosing several val¬ 
ues of e and z_ in the range (0,1) and numerically solve 
the equation Wr/wg = 1/2, 2/3 or 3/4 for p, given a fixed 
spin value a. The root finding and function evaluations 
required to solve the implicit equation for p are com¬ 
putationally expensive, particularly when using the first 
direct numerical method. Rendering Fig. [5] in Mathe- 
matica took two days on an Intel i7 PC with 8 cores each 
with 2GB of RAM using the first method. The second 
semi-analytic method that takes advantage of the com- 
mensurability of the frequencies and is also implemented 
in Mathematica, yields the same result in just under an 
hour. (A further speed up is likely for a Matlab or C-I--I- 
code). The results from the two methods agree and serve 
as independent checks. 

The numerically computed resonant surfaces are plot¬ 
ted in Fig.[5l This figure shows several surfaces or sheets 
corresponding to different black hole spins. The sheets 
are plotted on the {e,p, z'^} coordinate axis to make the 
correspondence to Fig. [3] clear. The blue sheets (a)-(c) 
indicate the location of prograde orbits. The spin value 
is decreased from sheet (a) to (c). The red sheets indi- 
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FIG. 6. The numerically computed location of the 1/2, 2/3 and 3/4 resonances for fixed eccentricity, e = 0.6 for the spin 
values a = 0.2, 0.5, 0.99. (These are the same spin values used in Fig. [5]) The resonances for prograde orbits are shown in 
solid lines and those for retrograde orbits in dashed lines. The resonances appear centered around the corresponding p* values 
given in Table H] with the 1/2-resonance (cyan) centered around p* = 8, the 2/3 resonance (blue and red) centered around 
p* — 10.8 and the 3/4 resonances (magenta) around p* = 13.7. The features of these plots agree qualitatively with those of 
Fig.E] 


cate retrograde orbits with spin values increasing from 
(d) to (/). The higher the eccentricity, the further out 
in p, the resonances occur. The more inclined the or¬ 
bital plane, the further out (closer in) the resonances 
for prograde (retrograde) orbits occur. These qualita¬ 
tive features are the same as those identified in the weak 
field limit in Sec. IV Bl Since the value of the orbital ec¬ 
centricity has the smallest impact on the location of the 
resonance, we subsequently choose a representative ec¬ 
centricity and explore the spin dependence for several 
low order resonances. In Fig. |6] the 1/2, 2/3 and 3/4 
resonances are shown for prograde and retrograde orbits 
with e = 0.6 around a Kerr black hole for several dif¬ 
ferent spin values. We see that the typical shape of the 
resonances shown in Fig. [5] is preserved. Prograde reso¬ 
nances with larger spin values are closer to the black hole 
in comparison to the point p* whereas retrograde reso¬ 
nances associated with higher spin values occur further 
outward. The shapes of these resonances are in qual¬ 
itative agreement with the weak held 'U'-'I' transition 
discussed with respect to Fig. |31 

A quantitative comparison between the numerically 
computed position for the 2/3, 6/7, 14/15 and 29/30 
resonance and the weak held approximation developed 
in Sec. IVBI is given in Fig. [71 The analytic result is ob¬ 
tained from Eq. (|32|), with e = 0.6 and varying spin val¬ 
ues. The weak held approximation provides a good ht for 
high orbital inclination (z^ = 1) even when p < 20. The 
weak held approximation deviates more strongly from 
the analytic results as the spin value increases, such that 
the orbits with a = 0.2 exhibit a better ht than those 
for which a = 0.99. For a given spin value the approx¬ 
imation hts the retrograde orbits (higher p) better than 


the corresponding prograde orbits (lower p). This is to 
be expected since in the approximation of Sec. IV Bl we 
used a quadratic approximation to the resonance condi¬ 
tion valid only for p —>■ oo and included spin effects only 
up to 0{a?/p^). Fig. [7] also shows the analytic solution 
on the equatorial plane as calculated from the roots of 
Eg. (1411) with the added second order correction for non¬ 
zero eccentricity p*e^/4(p* — 6) given in Eq. (14^ . These 
equatorial solutions are in very good agreement with the 
numerical results. They provide an easily computed in¬ 
dicator of where the weak held approximation strongly 
departs from the analytically computed surfaces. 


VI. THE LOCATION OF RESONANT 

SURFACES TERMS OF E, PARAMETERS 

When conducting numerical investigations of orbits in 
spacetimes that are more general than Kerr it is often 
useful to plot a Poincare map at hxed E and Lz as a 
diagnostic tool to explore the break down of KAM-tori. 
This method of surveying the parameter space is not as 
well suited for describing the physical location of the or¬ 
bit, but it unambiguously generalizes to metrics where a 
Carter constant Q does not exist. In this section we de¬ 
velop some intuition for the features of resonant surfaces 
in E, Lz space. This picture will be invaluable as we dis¬ 
cuss departures from integrability and torus breakdown 
in the next section. 

In Fig. [8] we show the 2/3 resonant surface for a maxi¬ 
mally spinning black-hole, a = 1, computed using the sec¬ 
ond numerical method of Sec. IV El In this hgure the dots 
correspond to computed points on the resonance surface. 
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FIG. 7. Quantitative comparison between the numerically calcu¬ 
lated resonance surfaces and weak field approximation of Sec. IV Bl 
for the 2/3, 6/7, 14/15 and 29/30 resonances. The eccentricity is 
fixed at e = 0.6. The weak field approximation is computed using 
Eq. II32II and is shown for spin values a = 0.2 in red (inner arch), 
a = 0.5 in blue (middle arch) and a = 0.99 in magenta (the outer 
arch). The analytic results on the equatorial plane are shown in 
matching colored dots, with the numerical solutions overlayed using 
black dashing. 



FIG. 8. The 2/3 resonant surface for a maximally spin¬ 
ning black hole a — 1 represented in terms of (p, e, zi) co¬ 
ordinates. Horizontal lines indicate longitudinal values of 
Z- — {0, sin-I, l/v/S, cos 1}, the bent arches indicated ec¬ 
centricities of e = {0, i, i, |, 1} 



FIG. 9. The 2/3 resonant surface for a = 1 projected onto 
E, Lz coordinates 0 ]. Lines indicate constant values of 
z?. = {0, sin l/v/2, cos 1}, and e = {0, L, |, |, 1} corre¬ 
sponding to those in Fig. [S] The large asymmetry across 
the Lz = 0 line is indicative of the large spin value under 
consideration. 


The lines indicate fixed values of eccentricity and longitu¬ 
dinal angles at which the sampling took place: horizontal 
lines correspond to = {0, sin(7r/8), I/-\/2, cos(7r/8), I}, 
the bent arches correspond to e = {0,1/4, 1/2, 3/4, 1}. 
For each point (p, e, z^) on the resonant surface, the 
corresponding E, and values were computed using 
Eqs. (fTOjl and (fTT|l . Care was taken to compute the ap¬ 
propriate Wx root associated with the retrograde p_ or 
prograde p+ branch of the resonant surface under consid¬ 
eration. The outcome is shown in Fig. [9l Note the large 
asymmetry about the Lz = 0 line that is due to the fact 
that the black hole is maximally spinning. The low en¬ 
ergy spike that coincides with prograde orbits close to the 
equatorial plane with low eccentricity is easy to miss nu¬ 
merically when exploring the parameter space in terms of 
E, Lz parameters. This spike constitutes a set of param¬ 


eter values that potentially have important astrophysical 
implications: Most particles in thin astrophysical discs 
around rapidly spinning black holes are expected to be 
prograde and on nearly circular equatorial orbits. As the 
spin of the black-hole decreases the resonance footprint 
on the E, Lz plane will become increasingly symmetrical 
as is shown in Figs. [TUI and ITT] 

In Fig. [T^the E, Lz projection of the 3/4 resonance 
surface is given for a black-hole with spin a = 0.5. The 
higher order resonance surface has features that are qual¬ 
itatively similar to that of the 2/3 resonance of equal spin 
shown in Fig. HOI For the higher order resonances how¬ 
ever the E values are slightly higher and the Lz values 
span a marginally larger range. 
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Lz 


FIG. 10. The 2/3 resonance for a = 0.5 in E, Lz coordinates 
for z?. = {0,sin|,l/V^, cosf,l}, and e = {0, 1, 1, |, 1}. 
Note the decrease asymmetry across Lz = 0 line when com¬ 
pared to Fig. [9] which is maximally spinning. 



Lz 

FIG. 11. The 2/3 resonance for a = 0.1 in E, Lz coordinates 
for z?. = {0, sin 1/^2, cos |,1}, and e = {0, 1, 1}. 

Since the spin is almost vanishing the figure is much more 
symmetrical about the Lz axis than Fig llOl andlOl 



Lz 

FIG. 12. The 3/4 resonant surface for a = 0.5 projected 
onto E, Lz coordinates, for z?. = {0, sin l/v^, cos-1,1}, 
and e = {0, j, j,l}. The main difference between the pa¬ 
rameter space covered by the 2/3 and 3/4 resonance in Fia llOl 
is that slightly lower energies are reached for the lower order 
resonance and for the higher order resonance the Lz values 
have a marginally larger span, the other features remain qual¬ 
itatively similar. 


From the figures in this section it should be clear that 
the representation of orbital parameters in terms of the 
(p, e, z^) coordinates is better adapted to the geometry of 
the resonant orbit, making it generally easier to quantify 
the resonance surfaces using these variables. Once ob¬ 
tained the resonant surfaces projected onto the E, and 
Lz orbital variables can potentially be used to aid nu¬ 
merical exploration into Kerr-like spacetimes discussed 
in the next section. 


VII. ROTATION CURVES AND THE 
BREAKDOWN OF KAM TORI 

A number of groups have numerically explored the 
breakdown of integrability in stationary axisymmetric 
vacuum metrics such as the Manko-Novikov metric that 
reduces to the Kerr metric for a certain choice of param¬ 
eters p^l - fisl . [TtI im . One of the tools used to explore the 
non-Kerr spacetimes numerically is to plot the Poincare 
map of the orbital motion for a fixed choice of energy (E) 
and angular momentum {Lz). An example of a Poincare 
map exhibiting broken tori is given in Fig. 1131 Each 
closed curve in the Poincare map corresponds to an orbit 
that was started on the equatorial plane and given an ini¬ 
tial momentum out of the plane. Each intersection of the 
orbit with the equatorial plane generates a point on the 
Poincare map. The radial coordinate (here p) and corre¬ 
sponding radial momentum {pp) of these piercing points 
are plotted. The closed curves in this Poincare map in¬ 
dicate that for most initial conditions geodesic motion 
in the Manko-Novikov metric remains integrable or reg¬ 
ular. However, among the closed curves, a Birkhoff chain 
of islands can be seen (enlarged in inset), which indicates 
that for certain initial conditions the regularity of orbits 
break down. This breaking of regular structure is as¬ 
sociated with resonances in the fundamental frequencies 
describing the orbit as predicted by the KAM theorem. 

The second diagnostic tool that gives insight during 
a numerical exploration is the rotation curve. This is a 
plot of the LOrjoje frequency as a function of initial con¬ 
ditions for the Hamiltonian potential being studied. For 
an integrable system the rotation curve is a smooth func¬ 
tion of initial conditions. In systems for which the KAM 
tori have broken, there are plateaus in the rotation curve 
where the numerically computed rotation number 
remains roughly constant. This situation is depicted in 
Fig.m 

In the numerical investigation of the Manko-Novikov 
metric, the breaking of low-order resonant tori is ob¬ 
served most dramatically at the 2/3 resonance. The ap¬ 
parent dominance of the 2/3 resonance over other integer 
ratios could heuristically be explained by the fact that 
in these studies, the deviation from Kerr was mainly a 
quadrupole perturbation and thus roughly proportional 
to cos 20. This geometric dependence of the perturba¬ 
tion indicates that any 2/m resonance is expected to be 
strongly excited. Integers other than n = 2 would arise 
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FIG. 13. A Poincare map generated for the Manko-Novikov 
metric. The inlay shows a close-np of one of the islands in the 
Birkhoff chain of multiplicity three. This island is associated 
with the breakdown of the 2/3 resonance as indicated using 
the rotation curve in Fig |14I The spin value for this simulation 
was chosen to be a = 0.9 and the dimensionless quadrupole 
deviation parameter q = 0.95 was chosen to be very large. 
The other orbital parameters are E = 0.95, Lz = S and rest 
mass p = 1. 



FIG. 14. The rotation curve is obtained by calculating the rotation 
number of each orbit from the Poincare map given in Figure 1131 
The rotation number, uq, which is equal to the ratio of orbital 
frequencies, is plotted as a function of the initial radial 

coordinate of the orbit. The rotation curve is generally smooth and 
monotonically increasing. However, corresponding to the Birkhoff 
chains of islands seen in the Poincare map, there is a plateau in the 
rotation curve for which = 2/3 is resonant. 


from the nonlinear coupling between the angle variables 
and the coordinates, which is a higher order effect. Since 
ujr < ujg for Kerr we always have that to > n, so 2/3 is 
the resonance with the lowest order integer ratio. It is 
thus expected to dominate the breakdown of integrability 
for a system subjected to a quadrupole perturbation. 

The numerical explorations of the breakdown of KAM 
tori has to date been limited to a very small subset of 
the allowed parameter space. To guide future studies 
covering the entire parameter space associated with the 
2/3 and 3/4 resonances discussed in Sec. IVII we now 
describe how to use the machinery developed in this pa¬ 
per to analytically compute the rotation curve for the 
Kerr metric. An important caveat to bear in mind when 


Rotation Number 



FIG. 15. The rotation curves for geodesics with E — 0.95^, in 
a Kerr black hole spacetime with spin a = 9/10. The differ¬ 
ent curves correspond to different Lz values uniformly spaced 
between Lz = 2.482p (lower curve) and Lz = 3.343/i (upper 
curve). For reference purposes, the 2/3 resonance value is in¬ 
dicated by the central red-line, the 1/2 resonance by the lower 
black-line and the 3/4 resonance by the upper black-line. The 
radius (r) refers to the radius of closest approach of the orbit 
to the black hole 


exploring large deviations from the Kerr metric is that 
although the torus structure is not destroyed it may be 
distorted and the breakdown of quasi-periodic orbits may 
be shifted from the exact location of the resonance pre¬ 
dicted by the rotation curve for orbits in Kerr spacetime. 
To zeroth order, however, the Kerr rotation curve pro¬ 
vides an indication of the interesting regions in the pa¬ 
rameter space. 

For the Kerr metric the rotation curve is computed by 
evaluating Eq. (ISU)) for a fixed E and using the pa¬ 
rameters defined in Eq. (1^51) . A representative example 
is given in Fig. [151 For a given E and L^, bound orbital 
motion is usually restricted to a small region of physical 
space. Equatorial orbits {Q = Z- = 0) have the largest 
radial extent while circular orbits (e = 0) have the largest 
longitudinal reach, z_max = z_(e = 0). These two ex¬ 
treme points define the boundaries of the rotation curve. 
We examine them more closely by first computing the 
variables {p,e,Z-} and subsequently the ratio ujrlujg. 

On the equatorial plane, Z-eq = 0 and the radial po¬ 
tential Eq. has one zero root, thus Wx = 0. The 
largest two roots of the remaining cubic determines Peq 
and Ceq. Eirst consider the effect of eccentricity Ceq where 
the two extreme cases for bound orbits are Ceq = 0 and 
Ceq = 1. Making use of Eq. (HDl) we have that if Ceq = 0 
the rotation curve reduces to single dot at 


r/ujg = 


I {p — 6)p — ia? ± Sa^Jp 
p^ + 3a^ =F 4a^ 


(51) 


If the E and values have been chosen to lie on the 2/3 
resonance, this value would correspond to either the lower 
rightmost or lower leftmost points in Figs. lU |T0] and 
ITT] depending on the spin of the black-hole and whether 
the retrograde (left) or prograde (right) orbit is under 
consideration. Eq. m in conjunction with Eqs. (1361) 
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and (l38l) can be used to estimate the characteristic E and 
associated with the k resonance in the non-spinning 
limit to be 

^ ^ 2(2 + ^)^ M ^ 

/r2 9(1+ k)’ ^2 > 

For the 2/3 resonance this evaluates to £1 = 0.95879/r 
and Lj = ±3.86702/r. If the spin is maximal, a = 1, the 
prograde branch has E = 0.88009/i, = 2.26471/i, p = 

3.86704 while the retrograde branch has E = 0.971899/x, 
Lz = —4.6614^, p = 15.733. These values can be com¬ 
pared to the lower corner points in Fig. |9l 

For parabolic equatorial orbits, Ceq —>■ 1 and E = p. 
This case corresponds to the largest possible radial ex¬ 
tent a marginally bound orbit can have. For the 2/3 
resonance these values are indicated by the upper left¬ 
most and upper rightmost points on Figs. IHl [IHl and [TT] 
We define the dimensionless angular momentum param¬ 
eter lz = Lz/p and note that on the equatorial plane, 
Wx = 0, n7+ = /z/2. Using Eqs. (flUl) and (fTTl) implies 
that the values of p are restricted to 



the given E and Lz- Substituting these results into the 
remaining expansion parameters of Eq. (1251) and eval¬ 
uating Eq. (15011 yields the maximum value the rotation 
curve attains. 


VIII. ESTIMATING THE SIZE OE THE 
PERTURBATION THAT RESULTS IN A 
DRAMATIC CHANGE OF DYNAMICS DUE TO 
TORUS BREAK DOWN 

The KAM criterion for the possible destruction of tori 
can be augmented by further estimates. A generic ana¬ 
lytic perturbation F{r, 0,pr,pe,P4,) to geodesic quantities 
such as the energy can be written as a Fourier expansion 
of the form 

E{r,e,pr,pe,P 4 >) = ^ E„(e,p, (56) 

n= —oo 

where q are the angle variables corresponding to the fre¬ 
quencies via 

+ (57) 


The factor under the square root is positive only if 
lz > 2(1 -I- VI — a) or lz < -2(1-1-Vl + «)• The smallest 
value of p that can be attained for equatorial parabolic 
orbits is thus = 2(1 -|- Vl T for prograde (p+) and 
retrograde (p_) orbits respectively. 

The parameters that enter into Eq. (E^ for equatorial 
parabolic orbits are 

yi 3p-ll ^1 ll-P ^2 ^ 

02 Vl li - 3p y2 


Evaluating Eq. dsni) we find that is a monotonically 
increasing function of p or equivalently lz- In the limit 
of a spin zero black hole the resonance = 2/3 occurs 
when p = 11.3273, lz = 4.18461. Parabolic orbits with 
larger \lz\ will not sample the 2/3 resonance regardless 
of inclination. For a maximally spinning black hole, the 
2/3 resonance occurs at p_ = 16.2914, lz = —5.01845 for 
retrograde orbits and at = 4.14634, lz = 2.53177 for 
prograde orbits (see Fig. IH]). 

Finally we consider the maximal vertical extent an or¬ 
bit can have. This occurs for circular orbits given a fixed 
E and Lz- If e = 0, 5i/yi = 0 and w+ = -2p+ 

Eq. (fTTl) gives Wx to be 


Wx = 


“ 2^) -f ‘i.a^lz + P (^(P - 3)p - (p- 2)2^ 

(55) 


which can be substituted back into in Eq. to 

give a polynomial in p that can be solved to find Pe=o for 


For most orbits, all the Fourier components in Eq. (1561) 
with n V 0 will be oscillatory and their contribution ap¬ 
proximately averages out over an orbit. The orbit’s lead¬ 
ing order secular evolution is driven by the component 
with n = 0. However, at a resonance where k • a; = 0 the 
Fourier components in Eq. (I56|) that involve the resonant 
combination (k • q) momentarily cease to be oscillatory 
since their phase becomes stationary. These components 
thus generically contribute order unity corrections to the 
secular evolution over the resonance time, except when 
the amplitude of these components is 0 (e) or smaller. 
Since for analytic functions the Fourier amplitudes £k 
fall off exponentially with k, we expect the resonance to 
have an appreciable effect on the dynamics only when 

Ok = El^*I^O(|ln(e)|). (58) 

i 

This criterion on k in (I58|) for “essential” resonances [s^ 
refines the bound for sufficient irrationality discussed in 
Eq. ([2]) for the special case of analytic perturbations. 

For EMRIs with the only source of perturbation being 
the gravitational self-force, even the low-order resonances 
are expected to be weak in the sense that the dominant 
dissipation Eq is generically larger than the resonance 
potential related to Fk. The resonances will therefore 
appreciably modify but not destroy the object’s contin¬ 
ued in-spiral. Now, however, consider the case of addi¬ 
tional perturbations that lead to very strong resonance 
modifications, where heuristically Fk Fq. Chirikov’s 
resonance overlap criterion specifies the conditions 
under which the complete loss of quasi-periodic motion 
is expected to occur as follows. Each strong resonance 
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captures orbits in its vicinity and is thus surrounded by 
an oscillation region similar to the Birkhoff islands shown 
in Fig. (11311 and manifests as a plateau in the frequency 
evolution similar to that in Fig (IH)) . The onset of full¬ 
blown chaos occurs when these Birkhoff chains associated 
with nearby resonances become large enough that they 
touch. More precisely, Chirikov’s criterion states that the 
transition to stochastic behavior arises when two neigh¬ 
boring strong resonances overlap in the sense that the 
frequency width of their oscillation regions is larger than 
their spacing in frequency. An estimate for the width of 
the resonance regions is [H, 

A|k • a)| -- -/e^Kk • Fhk|, (59) 

where Ja = {e,p^Z-) and Fkk are the forces such that 
Ja = f-Fa + O(e^) and Fa is Fourier expanded as ex¬ 
plained above. Overlap occurs for two resonances associ¬ 
ated with lattice vectors k and k' respectively when 

(60) 

for each of the frequencies uJi evaluated at the resonances 
of the unperturbed system. This estimate of the criterion 
for overlap is only a crude indicator for the transition to 
stochasticity, and in cases of interest the local dynamics 
must be systematically analyzed [^, • 

For most systems explored numerically to date, the 
plateaus in the rotation number curve remain small, and 
thus do not satisfy Chirikov’s criterion. This indicates 
that we expect weak chaos in the sense that torus dis¬ 
ruption is limited to a small region in phase space. 

IX. ASTROPHYSICAL IMPLICATIONS 

For astrophysical processes near super-massive black 
holes the Kerr geometry usually provides the background 
stage on which several small perturbation effects play 
their part. In select regions of the spacetime correspond¬ 
ing to low-order resonant orbits the smooth distortion 
due to perturbation induced effects is disrupted. The 
background spacetime geometry largely sets the location 
of these resonance induced disruptions imprinted on the 
dynamical structures of the perturbed system, while their 
details depend on the properties of the perturbations. 

Generally these effects are expected to be small, al¬ 
though persistent and robust. They are induced by all 
non-integrable perturbations of the Kerr metric. If the 
perturbation is not too large resonant effects occur at 
the locations quantified in this paper irrespective of the 
source of the perturbation. Table U indicates the charac¬ 
teristic length and timescales associated with resonances 
in dimensionless units and as well as in units charac¬ 
teristic of the region around Sgr A*. In Fig. [16] we 
have plotted all the low order resonances with order 
Ok = n m < 15 superimposed on an embedding di¬ 
agram for a non-spinning black hole. This plot gives an 



FIG. 16. The location of low order resonances around a black 
hole 0. Here low order resonances are plotted superimposed 
on a embedding diagram to give an indication of how strongly 
the spacetime in their vicinity is curved. The line-width in 
each case is scaled with 3 /Ok where Ok = m -\- n to give an 
indication of the relative importance of a particular radius. 
This image gives an over estimate of the importance of higher 
order resonances since the correct scaling according to Eq. m 
is K/Ol. 


indication of how strongly the spacetime for a particular 
resonant surface is curved and the relative extent of the 
regions of influence of the low-order resonances. Note 
the accumulation of resonances near the ISO which can 
also be observed in Fig. |4| and persists for the spinning 
case. In Fig. [T6| we have scaled the width of the lines 
demarcating each resonance by K{F)IOk with an arbi¬ 
trary choice of AT = 3 to result in good rendering. Recall 
that Arnold’s criterion, in Eq. (ED, governing the per¬ 
sistence of tori scales as KjO^. Higher order resonances 
are thus likely to be strongly suppressed compared to the 
schematic representation given here. 

The three lowest order resonant tori, 1/2, 1/3, 2/3 
whose parameters are indicated in bold in Table |T| are 
mathematically most likely to break. Of these it is 
likely that astrophysically the 2/3 resonance will have the 
greatest probability of being directly observable based on 
the following arguments: (i) For an EMRI, the impact of 
the broken resonance is likely to be larger if the perturb¬ 
ing object remains in the vicinity of the resonance for 
a considerable amount of time. The 1/2 and 1/3 reso¬ 
nances lie at 4i?5 and 3ARs respectively, close the ISO 
at 3i?s where the radiation reaction force is large and the 
orbit is transitioning to the plunge phase. The time and 
number of orbits that an in-spiralling object of mass /r 
spends near a given radius roughly scale as t ^ p^/{fi/M) 
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and N ~ j thus the influence of the 2/3 res¬ 

onance can accumulate loireer than for the 1/2 and 1/3 
resonances. See also Ref. [9| for more precise estimates. 

(ii) For electromagnetic observations, the fact that the 
1/2 and 1/3 resonances lie in the region of high curva¬ 
ture and possibly near or on the edge of the accretion 
disc is likely to confuse any distinct signal originating 
from this region. The 2/3 resonance on the other hand is 
located at roughly 5.4i?s further out of the potential well, 
making it easier for electromagnetic radiation to escape. 

(iii) The final reason for expecting the 2/3 resonance to 
dominate is akin to the discussion in Sec. IVIII where it 
was argued that quadrupolar gravitational potential per¬ 
turbations will preferentially excite the 2/3 rather than 
lower order resonances. 

If the resonance condition is satisfied the particle mo¬ 
tion either passes through the resonance, or, when the 
resonance dominates over dissipation so that Fk > Fq in 
Ea. (l56l) holds and depending on the initial conditions, the 
object can temporarily be captured in the resonance. For 
captured particles the resonance effectively fixes the or¬ 
bital frequency ratio of the object’s orbit which manifests 
as a plateau in its rotation curve, as seen in Fig.[T31 This 
is indicative of a stable, attracting resonance. Repelling 
resonances, typified by an inflection point or jump in the 
rotation curve, are associated with unstable periodic or¬ 
bits [1^, and will be short-lived in the presence of dis¬ 
sipative effects. The detailed study of the dynamics when 
a small mass enters the resonance region under various 
forms of perturbation is the subject of future work and 
has been studied in a particular Newtonian case in [^ . 
In the event of resonant capture the orbital parameters 
will linger within the resonance surface, possibly altering 
constants of motion because of the interchange of energy 
and angular momentum between the perturbation and 
the orbit. When the gravitational radiation reaction is 
significant, the dissipation will most likely cause the orbit 
to evolve towards a lower energy state, i.e. increasingly 
circular equatorial configurations, corresponding to the 
lower right-hand corner of Figs. |H]and|ni 

Now consider a collection of particles in a low-density 
accretion disk around an astrophysical black-hole. As 
explained in the introduction, their orbits will be influ¬ 
enced by a number of small perturbations which could 
cause them to become captured by a resonance, since 
the gravitational dissipation is small for low mass ratios. 
The evolution towards prograde circular equatorial orbits 
during resonance will limit the collisional interaction of 
the swarm of trapped particles, potentially creating a 
cohesive structure. It is thus possible that the charac¬ 
teristic structure associated with the resonances, such as 
the resonance zones visualized in Fig. 1161 will be im¬ 
printed on any thin disk surrounding a black-hole in the 
form of density inhomogeneities in much the same way 
as the resonant structures imprinted on Saturn’s rings 
Q. Unlike the rings of Saturn where matter largely re¬ 
mains captured indefinitely, the black-hole rings will be 
dynamical because radiation reaction will alter the res¬ 


onant structures and enable escape from resonance (as 
can be quantified using e.g. the methods in (6^). In 
this scenario, when a trapped over-density becomes mas¬ 
sive enough for the radiation reaction force to become 
important the black-hole ring will partially disintegrate, 
depositing an over-density of matter on the next inward 
ring in the disk. There, the matter will again be trapped 
for some time before continuing the in-fall. Thus in one 
related catastrophic event the whole ring structure will 
re-adjust, with the emitted radiation carrying an imprint 
of the particular resonances involved. Provided the accre¬ 
tion rate is slow enough the ring structure will reform af¬ 
ter each radiation reaction dominated ring collapse event. 
Since the radiation reaction force scales with mass ratio 
it is further expected that there will be a segregation 
in the particle sizes found in each ring. The tendency 
for captured orbits to evolve to a more circular, equato¬ 
rial configuration is expected to minimize the ejection of 
simultaneously captured particles due to collisional inter¬ 
actions. 

Another mechanism that could excite the resonant ring 
structure and lead to a distinct signature in the emitted 
radiation is the collision of a compact object with the 
matter in successive rings. The resulting collisional hot 
spot of excited gas will rotate with an azimuthal fre¬ 
quency set by the characteristic resonance surface. A 
detailed study of the possibility of a ring structure, its 
dynamics and observational signatures is left for future 
work. 

We now conduct a very precursory search for phenom¬ 
ena that could possibly be associated with the resonant 
structure around black holes. The closest super-massive 
black hole at the center of our galaxy presents an inter¬ 
esting test-bed. Recent monitoring of Sgr A* with the 
1.3mm VLBI showed time-variable structures on scales 
of ^ 4 i?s [^, . The physical origin of this structure 

is not yet clear, but it is interesting to note that the 
scale is similar to that of the low-order resonances given 
in Table I. As discussed above it is possible that tempo¬ 
rary capture of material or gas near the resonance loca¬ 
tion could lead to a time-varying signature in the photon 
emission. For argument’s sake, suppose that the origin 
of the structure at ^ 4i?s = 8M is in fact the 2/3 res¬ 
onance but displaced from the non-spinning value listed 
in Table [J because Sgr A* has spin. Using Eq. (H5|) the 
spin displacement of the prograde 2/3 resonant surface is 
p_i_ = 10.8 —5.36a. From the amount of spin displacement 
of the resonance needed to match the observed structure 
one could then infer that Sgr A* has a = 0.5. The plausi¬ 
bility of identifying this structure with the 2/3 resonance 
could be confirmed if the variability has characteristic 
timescales of slightly less than an hour. The increase in 
sensitivity of the VLBI measurements will enable resolv¬ 
ing more of the horizon scales, and it will be fascinating 
to see if the resonance structure can be revealed. Table U 
provides a quick reference for the possible characteristic 
time and length-scales. Note that because the coefficients 
in Eq. (1451) differ for different resonances, observing more 
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than one resonant structure places an independent check 
on the veracity of this method of determining the spin, 
since the displacement of both resonances due to spin 
must be consistent. 

Another example of a phenomenon that could poten¬ 
tially be associated with the orbital resonances is the 
quasi-periodic oscillations (QPOs) observed in the X-ray 
spectra of several black hole candidates. Four stellar 
mass black hole systems exhibit quasi-periodic variability 
with peaks at harmonic pairs of frequencies in a 2/3 ra¬ 
tio, one also shows an additional 3/5 ratio [i^. Recently, 
QPOs have also been identified in a super-massive black 
hole [i^ and in a tidal disruption event . A definitive 
physical explanation for the origin of the QPOs is cur¬ 
rently lacking. Numerous models have been proposed, 
including orbital resonances of any combination of the 
three orbital frequencies and their corresponding beats 
71,7%, accretion disk oscillations with nonlinear effects 
731 174 II. or variations in the geometry of the accretion 
flow [^. For the case where both the 2/3 and 3/5 ratios 
are observed, it is very interesting to note from Table H] 
that these two resonances are in fact nearest neighbors 
in p, with the 3/5 occurring at just slightly smaller p 
values than the 2/3 resonance. A disruption event at 
the 2/3 resonance could excite an event at the 3/5 res¬ 
onance because of their physical proximity. Measuring 
the frequencies of the observed resonances in QPO’s will 
give further clues as to whether they can correctly be 
attributed to orbital resonances around Kerr or whether 
other physics dominate over the orbital dynamics. 

The concentration of low order resonances near the 
black hole and their absence further out has important 
implications for the key science objective of testing the 
no-hair theorems using a super-massive black hole such 
Sgr A*. The no-hair theorems state that, provided the 
cosmic censorship and causality axioms hold, if the black 
hole’s mass and spin are known the quadrupole moment 
is fixed. Liu et al. [t^ have shown that recording the 
time of arrival signals of a pulsar orbiting Sgr A*, with 
orbital period ~ 4 months, for several years with the 
Square Kilometer array (SKA) will enable us to measure 
the mass of Sgr A* to a precision of 10“®, the spin 10“^ 
and the quadrupole moment to 10“^. This will allow 
a definitive test of the no-hair theorems. The detection 
of a pulsar even closer to the central object could allow 
the extraction of additional multipole moments through 
long term monitoring, thus mapping out more and more 
details of the structure of the central black hole. 

The analysis in this paper assures us that orbits around 
Sgr A* with a period of the order of weeks to months are 
sufficiently far from the low-order resonances that the 
KAM theorem guarantees quasi-periodic motion and the 
persistence of invariant tori under perturbations. This 
result implies that frequency drifts computed using per - 
turbative methods based on averaging, as done in [77|, 
accurately describe the physical system that is effectively 
free of stochastic motion. It also ensures that tracking the 
regular motion of a pulsar will reflect a map of the grav¬ 


itational potential it samples. We conclude this section 
by mentioning in passing the relevance to future gravita¬ 
tional wave detectors such as eLISA and their potential 
to directly probe resonant dynamics. A detector sensitive 
to the frequency band ^ 10“^ — 10“^ Hz observing an 
EMRI near SgrA* MsgrA* ~ 4 x 10 ®Mq [t^, can probe 
mean radial distances ranging from the event horizon to 
- 50M 0. This overlaps with the region where the 
strongest resonances occur (see Table |T] ). The potential 
direct detection of gravitational wave emission from a res¬ 
onance transit is an exciting possibility. It does however 
underscore the necessity to carefully model and incorpo¬ 
rate resonant effects in the search templates. The loss of 
phase coherence as the small mass object passes through 
a resonance could potentially make parameter estimation 
difficult. 


X. SUMMARY AND DISCUSSION 

In this paper we have computed the location and asso¬ 
ciated timescales of resonant orbits in the Kerr spacetime 
and comment on the observational and mathematical im¬ 
plications. We have considered resonances between the 
two fundamental frequencies corresponding to the radial 
and longitudinal motion, leaving the study of the effect of 
the rotational azimuthal frequency to future work. Our 
results provide a complete survey of the parameter space 
of resonant orbits, together with simple expressions for 
locating the resonances valid in the strong field region. If 
resonant phenomena are observed these expressions could 
provide and easy way of determining the spin of the cen¬ 
tral black hole. We considered several examples of elec¬ 
tromagnetic observations at different wavelengths that 
could be related to the capture and escape of material 
from resonances. 

We have computed the resonant surfaces both in terms 
of generalized Keplerian variables (related to the orbital 
geometry) as well as projected onto the E and pa¬ 
rameter space (related to the spacetime symmetries) to 
help identify promising choices of parameters for numer¬ 
ical investigations of torus breakdown in resonance re¬ 
gions. We have further found an analytic expression for 
the rotation curve associated with the Kerr metric as a 
function of E and that can be used for comparing 
to numerically computed rotation curves associated with 
Poincare maps. We expect low order resonances such as 
the 2/3 and 1/2 resonances to have the strongest observ¬ 
able effects on orbital motion, electromagnetic emission 
and the phase of emitted gravitational waves. These res¬ 
onances occur in the strong field regions of the spacetime 
at a distance of < 5.4 from the black hole and are 
fairly widely spaced by ~ 1.4 i?s in the limiting case of a 
non-spinning black hole. 

According to the KAM theorem low order resonances 
indicate where in a dynamical system the transition to 
chaos is likely to occur first. Such a transition requires a 
sufficiently large perturbation that could arise from vari- 
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ous sources such as the non-Kerr nature of the spacetime, 
the internal structure of the probe, alternative theories of 
gravity or the presence of other bodies or gas. We sum¬ 
marized the KAM estimate and additional arguments by 
Arnold and Chirikov to assess which resonant tori are 
expected to survive under the perturbation and where 
dramatic changes in the dynamics could occur. Com¬ 
bining these general bounds with our results for the reso¬ 
nances indicates that there is a large region for mean radii 
50 Rs ^ P ^ 1000 Rs where resonance effects are negligi¬ 
ble but where the spacetime curvature is sufficiently high 
that multipoles of the central object have an observable 
effect on the motion. This will enable tests of the no-hair 
theorem unimpeded by drastic changes in the dynamics. 

More stringent bounds than discussed in this paper 
on the occurrence of strong resonances and the transi¬ 
tion to stochastic behavior require detailed studies of the 
dynamics near resonant tori for the different kinds of per¬ 
turbations. We noted in Sec. IVlllI that for EMRls the 
gravitational radiation reaction force, which is important 
in the region where low-order resonances occur, sets an 
approximate scale for the strength of the perturbation 
required to destroy the tori. In future work we intend to 
study the details of the breakup of the tori and quantify 
the strength of the perturbation required for observable 
consequences in each case. 

Regions where lower order resonances occur also ear¬ 
mark the most likely positions in phase space where av¬ 
eraging methods must be modified to account for the res¬ 
onances. If an orbit passes through a resonance it effec¬ 
tively acquires a sudden change in the frequencies whose 
magnitude depends on the initial conditions. For suf¬ 
ficiently strong resonances, orbits can enter a resonance 
region, linger near it and subsequently escape. This man¬ 
ifests as a plateau in the frequency evolution similar to 
that seen in the rotation curve of Fig. [14] but tilted be¬ 
cause of radiation reaction. More work will be needed to 
determine how well the frequencies must be resolved to 
detect such a plateau and exploit the measurements to 
determine the system parameters with LISA. 

In addition to exploring the parameter space of res¬ 
onances in this paper, we characterize in the Appendix 
the locations of the innermost stable orbit beyond which 
geodesics plunge into the black hole. This provides a use¬ 
ful benchmark for the resonance locations and also indi¬ 
cates the region where zoom-whirl behavior would occur 
and where higher order resonances accumulate. 


Appendix A: Location of the Innermost Stable Orbit 


In this appendix we explore the location of the inner¬ 
most stable orbit (ISO) using the notation and variables 
given in Sec. IIII Bl Observationally this is an important 
location in the black hole spacetime since it demarcates 
the radius from which the test mass will enter the plunge 
phase of its orbit. Near the ISO a non-circular orbit ex¬ 
hibits zoom-whirl behavior as the radial frequency goes 
to zero while ujg and remain finite. The test mass will 
linger near its periastron for many periods of the (j) and 
9 motion, then rapidly zoom out to apastron and back, 
giving rise to a characteristic signature in the emitted 
gravitational waves. When tracking the frequency evo¬ 
lution of a particular orbit, this determines the final set 
of uj 0 frequencies that will be recorded before the 
test mass begins its plunge into the black hole. When 
discussing resonant orbits of a spinning black hole it is 
often useful to view their location relative to the location 
of the ISO for the same spin, inclination and eccentricity. 
The ISO surface (where vanishes) shares many of the 
qualitative features of the resonance surfaces (where the 
combination of frequencies muir — nuJe vanishes) discussed 
in the main body of the text and both are related to de¬ 
generate structures in the phase space that are broken 
under small perturbations. 

The ISO occurs when the potential barrier in the radial 
potential no longer shields the orbit from the singularity. 
Mathematically the condition that defines this orbit is 
that the middle two roots of Fq.(|3]) coincide. 

This allows us to express the and variables as. 


G7+ = r4 



Wx 


rj p 
1 -I- e' 


(A2) 


Fq. (1141) allows us to find an expression for in terms 
of r 4 , 


.2 ^ p(r 4 (e- 3) - {e + 3)p) 

2a^(e — l)(e -I- 1)^ 

p^(e — 3)^ r 4 ^ -I- 2 (e^ -I- 7 )pf 4 -I- (e -I- 3)^p^ 

^ 2a2(e-l)(e + l)2 ’ 

(A3) 
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where we have set rq = (1 -|- e)r 4 to yield a more compact 
expression. Eq. (fT^ provides a quadratic equation for 
r 4 that allows us to determine its value in terms of p and 
e. The solution can be written down in the form, 

1^4 = { — Bri + \/^j^)/4lr4, (A4) 

with 
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FIG. 17. Graphical representation of ISO surface constructed using Eq ’s. (IX3l)(IMl) . The ISO surfaces displayed in the top two 
images shows the same qualitative 'U' V' I' transitions as the resonant orbits. Smaller images in the bottom row give 
the top view of the resonant surfaces for various spins and eccentricities. The top left and the two images on the bottom left 
display the eccentricity dependence for different spin values. The ISO is linearly dependent on eccentricity (as opposed to the 
quadratic dependence observed in resonant surfaces). The spin dependence at various fixed eccentricities is illustrated in the 
right three plots where the typical 'V' profile is obvious irrespective of eccentricity. 


Ar^ = 64(e + 1)/ [p-{e- 1) (d - 1)] [p + {e - 1) (d + 1)] [p - {e + 1) (d + 1)]^ [p+{e + 1) (d - 1)]^ , 

Ar^ = a'^{e - l)^(e + 1)^ - 4a^(l - e)(3 - e)(e + l)^p + 2(e + (a^(3 + e)(e - 1) + 2(e - 3)^) 

- 4(e^ + 7)(e + l)p^ + (e + 3)^^"^, 

Br^ = a^(e — 3)(e — l)(e + l)^p — 4a^(e + 1)^ (e^ — 2e + 3) + 2(e + l)p^ ~ 5) + 2e^ + 14) 

- 4 (e^ + 2e + 3) + (e + 3)p®. (A5) 


We have set d = Vl — a^. In Eq. (IA4I) we chose the ‘+’ 
root of the quadratic to yield the correct behavior for . 
The ISO surface exhibits many of the characteristics seen 
in the resonance surfaces studied in Sec. IS These include 
the fact that as inclination increases the p values of the 
prograde/retrograde ISO’s respectively increase/decrease 
until the two branches coincide for = 1. This behav¬ 
ior is due to the fact that polar orbits are less influenced 
by the spin of the black hole. Also note that for = 1, 
dz^/dp = 0. A graphical representation of the ISO sur¬ 
face is given in Fig. [T71 


We now give a number of easily evaluated formulae 
for special parameter values. For circular orbits around 
a maximally spinning black hole (e = 0, a = 1), the 


behavior of as a function of p is described by, 

2 ^ P^ {-^Vp + 2a/3p -b 2^ + 3-3) 

3yp-l • ^ ^ 

Furthermore when e = 1, a = 1, 

2 SV2py^ - 3p^ - 4p3 + 4p^ 

4(9p2 + 4p + 4) ■ ^ ’ 

On the equatorial plane z_ = 0 and r 4 = 0. Just as 
for the resonances the strongest spin dependent effects 
can be observed here. The behavior of the ISO can be 
well characterized by merely examining the behavior on 
the equator. As inclination increases, the location on the 
ISO will lie between the extremes of the prograde and 
retrograde values found on the equatorial plane. The 
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characteristic 'V- 'V transition seen in the resonances 
occurs for the ISO as well. If r4 = 0 only the constant 
term (with respect to r4) in Eq. (fT^ with Eq. (IA2I) sub¬ 
stituted remains. This term results in a quartic equation 
for p that implicitly defines the location of the ISO on 
the equator, as also found in (t^ . 

(p — 6 — 2e)^p^ -I- 0^(1 -I- e)^(3 — ef 
-20^(1 -f e)p [2(6^ -f 7) -h (3 - e)p] = 0. (A8) 

If the black hole is non-spinning, a = 0, then p = 6 -I- 2e. 
If it is maximally spinning, a = 1, there are 3 roots to 
Eq. (IA8I) . namely p = 5 -I- e ± Ay/l + e and a double 
root at p = 1 -I- e. Of these the correct limiting cases 
for the retro and prograde ISO orbits are p_ = 5 -|- e -I- 
4VI + e and p+ = 1 -|- e respectively. Thus the maximum 
splitting between the pro and retrograde branches on the 
equatorial plane is p_ — p+ = 4(1 -|- ^/T+~e). 

If a 7^ 0 the ISO can be found by solving the quartic. 
(For plotting purposes it is easier to consider Eq. (jA8l) to 
be a quadratic equation for and plot the square root 
of the result as a function of p.) For convenience, the 
leading order spin dependence is also given here. 


Pt 


= 6 -I- 2e ± 4a 



0{a^). 


(A9) 


Appendix B: Properties of Carlson’s integrals 


Carlson has provided us with a number 

of symmetric, rapidly convergent schemes to evaluate el¬ 
liptic integrals both numerically and analytically. This 
appendix summarizes the identities and properties asso¬ 
ciated with Carlson’s integrals used in the body of the 
paper. Carlson showed that any elliptic integrals of the 
form 

h = r 

Jy yj (ai -I- /3it)(a2 + P2t){a3 + f5^t){aA + Pit) 

(Bl) 

can be expressed as 

(B2) 


where 


t/„- = + Y,Y,XuXm)l{.x - y), 

Xi = (ai + Yi = (ai + 

(B3) 


and 




1 /■°°_ (U _ 

2 Jo + a){t + j3){t -\- 7) 


(B4) 


This symmetric representation greatly reduces the num¬ 
ber of parameters from 10 in Eq. m to the three argu¬ 
ments in Eq. (IB4I) . All integrals have the same boundary 
conditions and the arguments enter only in the denomi¬ 
nator of the integrad. An added advantage of using Carl¬ 
son’s integrals is that they obey a number of identities 
that make manipulation of the parameters that enter as 
arguments possible without necessarily evaluating the in¬ 
tegral. These identities include the duplication theorem, 

Rpia, l3,"f) = 2RF{a + X, P + + X), (B5) 

where A = (a/3)^/^ -|- -I- the fact that 

Carlson’s function is homogeneous of degree — 

RpiXa, A/3, Ay) = A-^/2i7F(a,/3,7), (B6) 


for any A and a number of special symmetric cases that 
can easily be evaluated, 

Rf{I3, 13,13) = A-f(0, /3, /3) = ■ (B7) 

If the first argument is a = 0 a restricted version of 
the duplication theorem also holds: for a; > 0, y > 0, 
2; = {x + y)/2, 

Rf{Q,x^ ,y'^) = RF{Q,xy,z^). (B8) 

For small deviations from the symmetric case it is pos¬ 
sible to construct a rapidly converging series. We now 
derive this series for i?F(0, y -f <5, y — i5). Start with the 
integral representation of Rp, Eq. (lB4l) in the special 
case, 

RF{0,y + S,y-6) = ^ f JL, , 

2 Jo ^/t[t + y + ^)(t + y-^) 

(B9) 


and make a Taylor series expansion of the integrad in 
terms of S. Integrate the result term by term to yield, 

, r ^ ^ 3 ( 5 ^ 105 ( 5 '' 

+ = Vi(i + Tv + T024t 

TT / 1155(56 2252255® 

^ 2 ^ Vl6384y6 4194304y® ^ ' 

(BIO) 


Since only quadratic terms in the parameter S/y appear 
this series converges very rapidly. To complete the dis¬ 
cussion on Carlson’s integrals and mainly for comparison 
with other work we now give the relationship between 
Carlson’s integrals and some of the elliptic functions, 

i?/(0,l - fc^l) = A:(fc), (Bll) 

Rf( 0, a, /3) = /3-^/^K(^l - a//3). (B12) 

Here K(k) is the complete elliptic integral of the first 
kind. 
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